DCPT-02/53 



hep-th/0208202| 



Integrable aspects of the scaling (/-state Potts models II: 

finite-size effects 



Dept. of Mathematical Sciences, University of Durham, Durham DH1 3LE, UK 
2 IFT/UNESP, Instituto de Fisica Teorica, 01405-900, Sao Paulo - SP, Brasil 



We continue our discussion of the g-state Potts models for q < 4, in the scaling 
regimes close to their critical and tricritical points. In a previous paper, the spec- 
trum and full S-matrix of the models on an infinite line were elucidated; here, we 
consider finite-size behaviour. TBA equations are proposed for all cases related to 
021 and </>i2 perturbations of unitary minimal models. These are subjected to a 
variety of checks in the ultraviolet and infrared limits, and compared with results 
from a recently-proposed nonlinear integral equation. A nonlinear integral equation 
is also used to study the flows from tricritical to critical models, over the full range 
of q. Our results should also be of relevance to the study of the off-critical dilute A 
models in regimes 1 and 2. 



1 e-mail: p.e.dorey@durham.ac.uk 
3 e-mail: roberto . tateo@durham. ac .uk 




Abstract 



1 Introduction 



As a second-order phase transition of a lattice model is approached, its correlation 
length diverges. In the so-called 'scaling region' near to the transition, a continuum 
limit can be taken, and the model can then be investigated using the techniques of 
continuum field theory. The g-state Potts models [jl], |2| illustrate this notion very nicely. 
For q < 4 the addition of vacancies allows them to be defined so as to have both critical 
and tricritical points ||, and thus to have two distinct scaling regions, each with its 
own associated continuum field theory^. These are known as the scaling Potts models, 
with the words critical or tricritical being added if there is a need to be more precise. 

This paper is the continuation of a companion paper 0], in which the exact S- 
matrices of the scaling Potts models were discussed. Our starting-point there was some 
work by Chim and Zamolodchikov ||, who, noting that the models should be integrable, 
proposed a set of S-matrix elements describing the scattering of elementary kink-like 
excitations. Using the bootstrap technique and a 'minimal' hypothesis governed by 
the presence or absence of Coleman-Thun || type explanations of S-matrix poles, we 
were able to close the bootstrap for all of the critical scaling Potts models, and for the 
tricritical scaling Potts models with 4 > q > 2. 

The scaling Potts models can be related to c/>2i and 4>i2 perturbations of c < 1 
conformal field theories, and in this context an alternative set of elementary S-matrix 
elements had previously been proposed by Smirnov 0, using a construction based 
on reductions of the Izergin-Korepin ajj S-matrix [||. (Interesting features of these 
scattering theories have recently been discussed in pi IIO.) While the relationship 



between the two approaches is now being clarified 11], in ||] we found it more convenient 
to work entirely within Chim and Zamolodchikov's framework - the fundamental S- 
matrix elements and the formal vacuum structure are then continuous functions of q, 
and the connection with the Fortuin-Kasteleyn ||] formulation of the lattice model and 
its symmetries is rather more direct. Smirnov's approach is perhaps more natural if 
one wishes to discuss perturbations of specific minimal models, but, as we shall review 
below, the identification of the scaling Potts models with such perturbations hides a 
number of subtleties. Nevertheless, it does underline that a study via finite-size effects 
might be worthwhile, and this topic forms the main theme of the present paper. Taking 
as partial input the mass spectrum and S-matrix elements found in [Q], we propose sets 
of thermodynamic Bethe ansatz (TBA) equations describing the ground-state energies 
of the critical and tricritical scaling Potts models, for the values of q for which the 
associated minimal model is unitary. These proposals are checked in a variety of ways, 
and we also study aspects of the finite-size behaviour of the models using the so-called 
non-linear integral equation technique. 

The plan of the paper is as follows. Section ^ discusses the conformal field theory 
descriptions of the critical and tricritical points. Section summarises the necessary S- 



* Strictly speaking the direction in which the critical point is approached must also be specified - here, 
apart from in section 3.2, we only consider the (first) thermal direction, which is related to changes of 
temperature. 
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matrix results obtained in || and Q . Section ^| contains a short review of previous work 
on the TBA for foi and eft 12 perturbations, and discusses some features of TBA systems 
in general. We also sketch some of the reasoning which led to our main conjectures. 
The conjectures themselves are outlined in sections ||| and H in the form of sets of rules 
for the construction of the eight new families of massive TBA equations for fai and 4>\2 
perturbations of the minimal models _M Pi p+i, related to the critical and tricritical scaling 
Potts models at the particular values yfq = 2 cos(^j) (p = 6, 7, . . .) and ^fq = 2 cos(^) 
(jj = 5, p = 7,8,...) respectively. The first four sets of equations for the perturbed 
critical models are given explicitly in section ||, while the story for the tricritical models 
is illustrated in section ^ by the set of equations for the Ai^fi + 4>\2 theory. This is 
related to the tricritical branch of the Potts model at ^fq = 2cos(^). 

In section |7| our TBA systems are subjected to a number of analytical and numerical 
tests, all of which they pass. Further verification is provided in section ||, where a variant 
of the non-linear integral equation of |T^] is proposed to describe the finite-size ground- 
state energy at arbitrary values of q and its solutions are compared numerically with 
those of the TBA equations. This section also shows how the nonlinear integral equation 
technique can be used to study the interpolating flows between the tricritical and critical 
models, by taking an equation first introduced in []l3| and tuning its parameters to 
suitable values. 

The full set of TBAs related to both 021 an d (fti2 perturbations of minimal unitary 
models, and the associated sets of functional relations, are given in four appendices. 

We end this introduction with a remark on the possible wider relevance of our 
results. In the following, we have concentrated on the integrable quantum field theories 
associated to the continuum limits of lattice Potts models near to their critical and 
tricritical points. A renormalised field theory contains information which is universal in 
nature, and its relevance is not restricted to any specific member of a universality class. 
Other lattice systems associated with the "Potts" universality classes are the dilute A 



models of 14, 15]. While these lattice models are only defined at discrete values of q, 
they have the advantage of being soluble not only at but also away from the critical 
point, even on the lattice. The link with Potts models comes via the identification of 
their scaling limits with the 4>i2 and fax perturbations of unitary minimal conformal 



models. The exact correspondence is [16 



Dilute A p _i , regime 1 ± «-► M p . , p +i + 4>2i *-> y/q = 2cos(-£j-) , critical branch; 



Dilute A p , regime 2 <-> M. P:P +x + 4>\2 <-> y / g = 2cos(^), tricritical branch. 

These are precisely the points at which we have been able to conjecture TBA descrip- 
tions of the continuum models. Past experience (see []l7| , 18] for examples directly 



relevant to the matter in hand) suggests a close link between the underlying mathe- 
matical structures of the lattice and continuum models, when both are integrable. We 
therefore expect that many of the results reported in this paper, such as the general 
forms of the TBAs and the Y-systems, will also play a role in the study of the off-critical 
dilute A models, at least in regimes 1 and 2. 
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2 The conformal field theory description of the critical 
and trier it ical points 

If a model is placed precisely at a second-order transition, its correlation length is infinite 
and it has no intrinsic length-scale. Its behaviour should therefore be described by a 
conformal field theory jOJ. Some key features of the conformal field theories relevant 
to the critical g-state Potts models with q < 4 were identified by Dotsenko [^] (see also 
[ffill). His work made use of previous predictions for certain critical exponents of the 



(/-state Potts models |22|, 23, p4[ , in particular the following formula, first proposed by 
den Nijs || : 

Vt ^ = 3T2X Wh6re \/9 = 2 sm(^A) and < A < 3/2 . (2.1) 

Here, y T is the renormalisation group eigenvalue for the energy operator e , given in 

terms of the scaling dimension x t of e by y T = 2 — x e . It is related to the specific 
heat exponent a as a = 2 — 2/yx- (In fact, den Nijs gave his formula in terms of 
(jl = tt/2 — ttX/3, with y/q = 2cos/n, but the parametrisation in terms of A will be more 
convenient below.) 

At the same time, conformal field theories with central charge c < 1 can be parametrised, 
at least partially, by a real number £ > such that 

For any (not necessarily rational) value of £, these theories admit so-called 'degenerate' 
primary fields, which are special in that they have null fields in their sets of descendants, 



causing their correlation functions to satisfy differential equations [19|. The possible 
(left) conformal dimensions these fields araj 

A,„(fl - m %-^( - 1 (r,. 6 N), (2.3) 

with a similar formula for the right conformal dimensions A r / s / . The spinless degenerate 
primaries have A = A so that their scaling dimensions are x r>s = 2A r]S ; these fields will 
be denoted 6 rs . 



In [19] and J20| |, the 2- and 3- state Potts models were identified with c < 1 conformal 
field theories with £ = 3 and 5 respectively. In both cases, the energy operator e was 
shown to correspond to the primary field <p2i ■ In [ 20 1 , Dotsenko conjectured that the 



same should hold for all q < 4, implying the general relation 



y T (0 = 2-2A 2il (£) = ^!^ (2.4) 



f In ||,||, the conformal dimensions are defined such that A r , s becomes A s r 
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and, comparing with (2.1), 



3 + 2A _ 3 (g-1) 

An immediate check on this hypothesis comes via the operator algebra 021021 ~ ^+031 > 
which predicts that the critical exponent y T (the second thermal exponent) should be 
given in terms of the scaling dimension 23^ as 

4 4(3 -2A) . , 

y =2-x 3 i = -- = — ) -f- . 2.6) 

y T2 <M ^ (3 + 2A y ' 



This matches the value calculated by Nienhuis |23j using a mapping onto a Coulomb 
gas. Note, since < A < 3/2 for the Potts models, £ lies in the range [l,oo] and 
-2 < c(0 < 1. 

So far so good; but one should beware that the value of c does not specify a conformal 
field theory uniquely. Consider the situation when £ is rational, and suppose that 

£ = -r— (2-7) 

p — p 

with p and p' coprime integers with p' > p. Then fl2.2| ) and (^j^) turn into the familiar 
formulae 

_ -. _ W_ ~ P) 2 a (p'r - ps) 2 - ()/ - p) 2 
pp' App' 

and the degenerate fields can be identified in pairs as 

4> rs = 4>p- r , P '-s ■ (2-9) 

Operator product expansions between fields in the subset 

{4> rs \ r=l...p-l,s = l...p'-l} , (2.10) 

and their descendants, then close amongst themselves. The conformal field theory 
containing just these fields is called the diagonal, or 'A' series, minimal model M pp i ■ 
Depending on the values of p and p' , there may be other consistent truncations to 
other field theories also containing finite numbers of degenerate primary fields. All are 
examples of rational conformal field theories, and the full set of possibilities for c < 1 
forms the famous ADE classification of |25|| . However, apart from the special values 
£ = 3 and 5 (corresponding to q = 2 and 3) none of these is the conformal field theory 
of a g-state Potts model. This follows from the fact that the relevant torus partition 



functions do not coincide [26|. 



At first sight, this might seem to contradict the claim of [25] to have found a classi- 
fication of modular invariant partition functions with c < 1. The Potts model partition 
functions are certainly modular invariant, so how can they escape this result? The ex- 
planation is that the proof in [25] made essential use of the requirement, first emphasised 
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in [27], that all characters should appear in the partition function with non-negative 
integer multiplicities. This must be true for local quantum field theories whose partition 
functions can be given as traces of powers of a transfer matrix, but it fails for the Potts 
models at general values of q, for which no such transfer matrix can be defined. 

Another way to see the special nature of the Potts conformal field theories is to 
consider the limit q — > 4, or £ — » do. To take this limit through minimal models we 
must send p,p' — ► do, keeping \p— p'\ finite. The scaling dimensions of the degenerate 
fields tend to x rs = (r — s) 2 /2, as is reasonable for a theory with c = 1. However, 
notice that x T)S = x r+ k )S +k for any k . Taking the limit £ — > oo through a sequence of 
minimal models therefore results in a theory in which each scaling dimension appears 
with an infinite degeneracy^. By contrast, the degeneracies of the scaling dimensions in 
the 4-state Potts model are finite. 

The discussion so far has concentrated on the critical Potts models. For the tricritical 
models, it turns out that universal quantities are still given by formulae such as ( |2,1|) 
and (|2.2| ), with the same formula ^fq = 2sin(|A) relating q to A or £ = (3+2A)/(3— 2A), 
but with A now required to lie in the range 3/2 < A < 3 Q. At fixed q this is achieved 
by sending A to 3 — A, or £ to —2 — £. Under this continuation, c(£) is mapped to 
c(£+l), and A2i(£) to Ai2(£+1). This is in accord with the observation of fl2l| , that 



the exponents for the tricritical models 22, E3] can be recovered by the same style of 



argument as followed for the critical models, if the energy operator is identified not with 
02i, but rather with 0i2- However, just as for the critical models, it is only when q is 
an integer that the partition function of a tricritical Potts model coincides with that of 
a minimal model. 

For both the critical and the tricritical models, the conformal dimension A of the 
energy operator e is 

^—h+w^x)- (2 - u> 

Note that A(0) = 1, while A(3) = 0: the range [0,3] for A is precisely the range for 
which the energy operator is both relevant and of positive conformal dimension. 

We have stressed the fact that the critical and tricritical Potts model partition 
functions do not generally coincide with those of minimal models, even at points where 
the central charges match, because a knowledge of the operator content as encoded 
in the partition function is crucial for the proper interpretation of finite-size effects. 
Imagine that the theory is defined on a cylinder of circumference R. If it is conformal, 
then the only scale is provided by the system size itself, and so E(R), the ground state 
energy, must be proportional to 1/R. The precise relation is [|30|, [31]] 

E{R) = -^c efr (2.12) 



*This limit needs some care. Here we have implicitly focused on a finite set of ('Kac'-like |2q| ) 
operators, taken the limit q — > 4, and only then allowed the number of operators to tend to infinity. 
Taking the limit in the other way gives a theory even less like the 4-state Potts model |29(] . 
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where c e g, the 'effective central charge', is equal to c— 12x m i n , with x min is the smallest 
scaling dimension in the model. This is the dimension of some field (f> vac which generates 
the vacuum, or ground, state for the conformal field theory on the cylinder. For a unitary 
theory, there are no negative-dimension fields and 4> vac = I, the identity operator. 
Thus x m i n = 0, and c e fr = c. However, if the model is non-unitary then negative- 
dimension operators are to be expected, and as a result c e g- is larger than c. Whether 
the Potts models should be considered as unitary or not is perhaps a moot point, but an 



examination of the partition functions in 26 ] shows that all the scaling dimensions there 



are non-negative, and hence, for all < q < 4, for both the critical and the tricritical 
models, we have c e g = c. 

We are interested in the scaling regions around the critical and tricritical points. In 
the continuum limit, these should be described by perturbed conformal field theories of 



the form [32]: 



)d l x (2.13) 

where Acft is the action at the critical or tricritical point, r measures the (scaled) 
deviation from the critical temperature, and the energy density operator e(x) can be 
identified with 4>2i(x) or cpi2(x) for the critical or tricritical points respectively. The 
dimensionful coupling r introduces an independent length scale m oc t 1 ^ 2 " 2 ^ , where 



A = A21 or A12 is the conformal dimension of e, given in terms of A by ( 2.11 ). The 
ground state energy E(R) can still be expressed as in ( |2.12| ), but now the effective 
central charge can depend on the dimensionless quantity r = mR, and we should also 
allow for the possibility of a 'bulk' term in E(R), proportional to the system size: 

E(R) = £( r )fl- JL c ( r ). (2.14) 

The 'scaling function' c(r) (for brevity, we shall omit the subscript 'eff ') encodes a great 
deal of information about the off-critical model, and will be the main topic of this paper. 



3 Mass spectrum and S-matrix data 

An important input to our conjectures is the infrared information provided by the 
mass spectrum and S-matrix of the model on an infinite line. To make this paper 
relatively self-contained, in this section we summarise the relevant data from Q. The 
full spectrum contains both kink and breather states, but since only the diagonal S- 
matrix elements (those involving at least one breather) enter directly into the TBA 



systems, only these will be quoted here. The S-matrix elements in § |3.2| , §3.3 and £ 3.4 
depend on the parameter A introduced in the last section. For rational values of A, our 
results will be equally applicable to perturbations of minimal models M pp ', modulo the 
issues of choice of vacuum state 4> VSLC discussed above. (For the unitary cases p' = p+1 
which form the main concern below, the ground states agree and these do not arise 
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anyway.) The relation between A and p and p' is 

A = p2 — | , (02i perturbations) (3.1 

V _ 3 

p 2 ' 



A = ^ — | , (0i2 perturbations) (3-2) 



3.1 Perturbed critical models: < A < 1 (one particle) 

In this regime there is only the fundamental kink K = K\. The physical-strip poles 
in the kink S-matrix elements at 6 = ^ and 6 = % are due to the fundamental kink 



3.2 Perturbed critical models: 1 < A < 3/2 (two particles) 

A pole at 6 = m{l—l/X) enters the physical strip as A passes 1, and is interpreted as 
neutral bound state B = B\ with mass rriB 1 = 2m^cos(^ — ^). This extra particle 
brings with it two new scattering amplitudes, Sb 1 k 1 and Sb 1 b 1 '■ 

Sb iKi = + + ( 3 - 3 ) 

sb iBi = mm-i], (3.4) 

where 



nva 1 



H = (a)(l - a) and (a) = Sm *^f + . 2 ( . (3.5) 

smh(| - *ft) 

Even though Sb 1 k 1 and Sb 1 b 1 contain further poles, for values of A in this region they 
can all be explained Q by invoking a variant of the Coleman-Thun (fj, ^] mechanism. 



3.3 Perturbed tricritical models: 3/2 < A < 2 (four particles) 

In addition to K\ and B\ , two extra particles, a kink K2 and a neutral breather B3, 
now enter the spectrum. Their masses are 

m K2 = cos(§ - j^) , m Ba = 2m Bl cos(^ - f ) . (3.6) 

No further particles are needed to explain the pole structure up to A = 2. The new 
diagonal S-matrix elements are as follows: 

s BlK2 = wmw+m-w 

Sb 3 b 3 = [|] 3 [i] 3 [|-i] 2 [| + i][|-I][f-|] 

s BzBl = [i + ^hh + m-^u-m-^] (3.7) 
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3.4 Perturbed tricritical models: 2 < A < 9/4 (six particles) 

For 2 < A < | , two more breathers, B2 and -B5, enter the spectrum. The scattering 
amplitudes for the new particles become increasingly complicated, and greater reliance 
is placed on the Coleman-Thun mechanism to explain the pole structure. The new 
S-matrix elements in this region are: 



Sb 2 b 1 = 


[1 


_ j_p _ j_ipip _ 11 

2AJL3 2AJL2AJL2A 3J 










Sb 2 b 2 = 


[§: 




-a] 2 








Sb 2 b 3 = 


n 


1 [ 6 ] [ a elie a] [a 2] Is 


a] 2 








Sb 2 k 1 = 


n 


|[|][| + a][| + a] 










Sb 2 k 2 = 


[§ 


1 l 2 [i + 1 l 2 [i + sirs _ 

2AJ L6 1 2AJ L6 1 2AJL2A 


-1] 








Sb^B! = 


[i: 


| 2 [|-a][|+a][a-|][a- 


|][i-a] 2 








Sb 5 b 2 = 


[f 


- ^] 2 [l-^] 2 [§ + ^][§ + 


J_l 3 [l _ J_- 

2AJ 2A. 


i 2 [^- 


|] 


P - 21 
L2A 3J 


Sb 5 b 3 = 


[i 

Le 


- 1 IT! + 3 1T3 _ J_1 2 [A _ 
2AJ L6 ' 2AJ 12 2AJ L2A 


|][|+2\] 3 ! 


■ 3 1 

-2A 6. 


1 3 [ 


s _ lira 

2A 2 J L 6 


Sb 5 b 5 = 


[I 


-!] 2 [I] 5 [!-i][f-|][i] 5 i 


i+A] 3 [|- 


Iff!] 


[| 


-a] 2 


S B 5 K! = 


[I 


- _Ll 2 [i _ _L1 2 [I + JILL _ 

2AJ L6 2AJ L6 ' 2AJ L2A 










Sb 5 k 2 = 


it: 


I 2 [|] 2 [|-a] 2 [|-a] 3 [| + a 


m-m- 









The complete mass spectrum up to A = 9/4 is: 



m Kl 


= m 






m Bl 


= 2m cos( 


71 
1 


- — ) 

2A / 


m K2 


= 2m cosl 


'tt 
> 3 


- — ) 

2A I 


m B3 


= 4m cos( 


71 
1 


~ 2T) COS (2T 


mB 2 


= 2m cos( 


TT 

2 


-5) 


m B5 


= 4m sin( 


5) 





(A>0) 

(A>1) 

(A>3/2) 

(A>3/2) 

(A>2) 

(A>2) 



4 The thermodynamic Bethe ansatz 

Our aim is to probe the models (|2.13j ) via their properties in finite spatial volumes. The 
main tool will be the thermodynamic Bethe ansatz (TBA) approach, a technique which 
allows the finite-size ground state energy of a theory to be obtained via the solution of 
a collection of coupled nonlinear integral equations. 

If the set of TBA equations is to be finite, then we would expect that the unperturbed 
theory should be associated with a minimal model in some way (as recalled above, for 
non-integer values of q the critical or tricritical g-state Potts model is never precisely 
the same as a minimal model, but the operator subalgebras generated by the energy 
operator agree, and so the ground-state energies based on the 'conformal' vacuum state 
/ should coincide.) The precise form of the equations will depend on the particular 
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model being perturbed, and only a limited number of cases have been studied to date. 
Before describing our new conjectures, we shall give a brief summary of the earlier work. 



4.1 Earlier work 

The situation is summarised in table |l|; the cited papers can be consulted for further 
explanations. Many of these results can, sometimes with hindsight, be labelled according 



to the l g o g' 7 scheme put forward in [42], where g and g' indicate a pair of diagrams 
of ADET type. Since the A\ incidence matrix is zero, the TBA for the A\ o A\ case is 
trivial, reflecting the fact that the field theory of + 4>\3, the thermally-perturbed 
Ising model, is free. Two other cases, labelled (G2) and (F4), also have a Lie-algebraic 
interpretation, though the fact that these algebras are not simply-laced places them 
outside the set of g o g' systems. However, their form is in line with a more general 
set of Lie algebraic TBA systems discussed in |^3j filf . As was remarked in Q, there 
is an intriguing coincidence between the sequence {Ai, A2, G2, D4, F4, Eq, Eq, E7, Eg} 
which arises naturally in connection with the Potts models, and the 'exceptional series' 
of Lie algebras making up the last line of the extended Freudenthal magic square, 
as discussed by Deligne, Cohen and de Man, Cvitanovic and others [45, [If], [17], 45 



In [47], Deligne and de Man identified an extra member of the series, namely the 
superalgebra OSp(l\2). From our point of view this corresponds to M7.10 + <p2i, where 
the unperturbed CFT is the OSp(l\2)W x OSp(l\2)^ /OSp(l\2)W coset conformal field 
theory [Q]. The relevant TBA is the T\ o D3 system, the n = 3 case of the third line of 
table @. 

However, it is important to realise that all of the infinite sets of systems listed in table 
Q describe the ground states of perturbations of non-unitary minimal models, generated 
by negative-dimension operators and not the identity. These TBAs are therefore not 
directly relevant for the g-state Potts models, although one might hope that the addition 
of a suitable chemical potential would allow them to describe the Potts vacuum state 
as well (see for example |5(], [51], |34|). We will not discuss this possibility any further 
here, but it would be an interesting avenue to explore. 

The observation that a TBA system related to G2 might describe the ^21 pertur- 



bation of .Mg^o (A = |) was made in [39]. This point lies in the region where the 



analysis of Q suggests two particle types in the model, and indeed that is precisely 
what the TBA system of [||] predicts. This case will be discussed along with the new 
TBA systems below, and for now we merely stress that, just as for almost all the above 
examples, this TBA system was first obtained as a conjecture, only subsequently being 
checked to describe the claimed model. 

A more deductive approach can be found in a paper by Bazhanov and Ellem pF 



who discussed the 4>2i perturbations of A4^^ (A = A) and (A = ^), taking 

Smirnov's RSOS description of the massive scattering theories Q as their starting- 
point. Subject to some mild assumptions, they derived sets of TBA equations for these 
two cases. Crucial was the fact [^] that for these models the RSOS restriction on the 
allowed vacuum states coincides with the restriction imposed on adjacent spins in the 
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Minimal 
model 



Perturbing 
field 



Ceff(O) 



TBA 

system 



■M-n,2n-l 

M 

2n+l,4n 
M 2n+l,4n-2 
M 3 ,4 
Ma,5 

■Mg.io 
(c=l) 

■M 4j5 

^2,5 
A^2,2n+1 



<P21 
021 
021 

021 (=013) 

021 

021 

021 

021=012 

012 

012 

012 

012 

012 

012 



An 



3 

4 

_9_ 

10 



3/i. 



(n-2)(4re-3) 

n(2n — 1) 
(2n-3)(4ra-l) 

2n(2n + l) 
4(n-l)(2n-7) 
' (2n-l)(2n + l) 

1 
2 

J_ 

10 

4 
5 

14 
15 



22 

55 

6 
7 

j_ 

10 

1 
2 

— 22 

5 

2(n-l)(6n-l) 
2n + l 



(n-2)(2n+3) 
n(2n— 1) 

4n 2 +2n-3 
2n(2n + l) 

4n 2 -4 
4n 2 - 1 

1 

2 

_7_ 

10 

4 

5 

14 
15 



52 
55 


7 

_7_ 

10 

1 

2 

2 

5 



■P12 
,2n-l 

,2[n±l] 
^11 

hi 
hi 
i 

hi 
hi 
hi 
hi 
h% 

hn 



T\ o A 2n -i 
T\ o T 2n -i 
TioD n 
A 1 oA 1 

A 2 oAi 

(G 2 ) 
D40A1 

(ft) 
E 6 oA 1 

E 7 oAi 

E 8 oA 1 

T X *Ay 



Table 1: TBA systems for minimal models perturbed by 02 1 an d 0i2- For the systems 
of gog' type, the mass in the magnonic system should be placed on the first node (that 
furthest from the tadpole or fork for T n or D n ). [x] is the integer part of x. 



hard hexagon model [53]. Unfortunately this equivalence holds only for a subset of the 
whole family of models related to (as mentioned in [)35|| , another example, to which 



we shall return in §|6| below, is the theory fA^fi + 0i2)- Thus, the extension of the 
approach of [35] to further models faces considerable obstacles. 

One other piece of work should be mentioned at this stage, even though its ultimate 
conclusions appear to run counter to the findings that will be reported below. In the 
course of a detailed study of character and polynomial identities associated with general 
02i perturbations, Berkovich, McCoy and Pearce [54] discussed possible TBA equations 
for all unitary minimal models A^ PlP +i, which reduced to the result of [^] for p = 4. 
They were only able to specify the general form of these equations, and certain ('kernel') 
functions necessary for a complete conjecture were left undetermined. Nevertheless, at 
least on a naive reading, the equations of [|54|] appear to entail just a single massive kink 
in the model for all values of p. This is at variance with the results of Q, which for 
02i perturbations predict that there should be a kink and a breather for all p > 5. In 
later sections we shall propose some new TBA systems which are consistent with this 
spectrum. There may be room to accommodate both our results and those of [54] - the 
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closing remarks of [55| indicate one possibility - but we suspect that the final resolution 



will reside in some modification to the considerations of [ 54 1 to bring their TBA systems 
into line with ours, at least insofar as they concern 021 perturbations of unitary minimal 
models|. 

For now we shall leave this question unresolved, and proceed with a discussion of 
the general form that we would expect the TBA equations to take, assuming for the 
time being that the results of § do indeed provide a reliable guide. 

4.2 The general structure of TBA systems 

When all scattering is diagonal, the integral equations of the TBA follow directly from 
the bulk S-matrix. The ground state energy of the model on a circle of circumference R 
is written in terms of 'dressed' single-particle energies e a {9) (pseudoenergies) [|38| . These 
pseudoenergies solve a system of non-linear integral equations of the following form: 

N 



e a {9) = Rm a cosh0 -J2®ac* L c (9) . (4.1) 



c=l 



Here L c {6) = ln(l + e £c ( )), * denotes the convolution 

f*9{0) = — j f(e')g(6-8')d6', (4.2) 

and the 2 — > 2 S-matrix elements S ac (9) influence the equations through the kernel 
functions <$> ac (9): 

$ac(0) = ~i^ In S ac (9). (4.3) 
The number of pseudoenergies coincides with the number N of particle types in the 



original scattering theory. From (4JJ, they have the large 6 asymptotics 



e a (0)|«-oo ~ Rm a cosh9 , (a = 1,2, . . . ,N ). (4.4) 

If off-diagonal scattering is also involved the complexity of the method increases sig- 
nificantly, as it is necessary to perform an extra 'diagonalisation' step before the final 
equations can be written down (see, for example, for further discussion of this 
point). But once the dust has settled, one generally finds that the parts of the TBA 
equations associated with diagonal scattering are as before, while the contribution of 
each kink multiplet K is split into two parts. The first is described by a single pseu- 
doenergy £k(9) with an asymptotic of the type Q4.4j ), with rriK the common mass of 
all kinks in the multiplet. In addition, diagonalisation results in the introduction of a 
(possibly infinite) number of auxiliary pseudoenergies. These behave as 

£a{9) | e^oo ~ constant , a = N + 1, N + 2, . . . , (4.5) 



§ In this context, it is worth noting that there is strong evidence that further novel sets of character 
and polynomial identities can be found using the TBA equations proposed in this paper. We would like 
to thank Ole Warnaar for performing an initial check of this. 
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and can be associated with fictitious 'particles' transporting zero energy and zero mo- 



mentum (note that m a cosh# in (4.4) is the single-particle energy on the infinite line). 
These new particles are often called 'magnons', and can be thought of as constructs in- 
troduced to get the counting of states right. The diagonalisation is by no means trivial, 
and in the following we will instead make a conjecture, using the above considerations 
as our guide and borrowing elements from some previously encountered TBA systems. 

4.3 The steps towards the main conjectures 

Our TBA conjectures were found via a link with the TBA equations for perturbed 
Z]\r systems of ]56[| . Originally this emerged from a study of four previously- known 
012-related cases. These were: M.^,a + 4>i2 (related to E$), + <pu (related to 

£7), + 0i2 (related to Eq) and «Mio,ii + <j>i2 (related to F4). No hint of a common 
structure came by just looking at the first three models, but a consideration of ^6,7+012 
and 10,11 + </>i2 alone revealed some striking analogies with systems related to . To 
see this, let us compare them with the sine-Gordon TBAs at ^ = -g and §7 = 7- The 
four systems are depicted in figures ^ g, [3| and ||. (Figure ^ is the -E^-related system: the 
distinction between solitons and antisolitons is special to this, the case of the tricritical 
Potts model, and so the Eq diagram symmetry has been used to superimpose each pair 
of soliton-antisoliton nodes, resulting in the diagram shown. Likewise, figure Q is the D4- 
related point of the sine-Gordon model, and the D4 diagram symmetry has been used to 
superimpose the soliton and the antisoliton.) There is a formal similarity between the 
pairs of diagrams in figures [l] and [2] and in figures |3| and ||, the main differences being 
related to the fact that there is just a single soliton-antisoliton pair in the sine-Gordon 
model, and correspondingly just a single massive 'soliton' node in the -related TBA 
systems. 

The sine-Gordon TBA systems of figures and I are the first members of a series 
of models at ^ = l 2 , represented in figure [|, which have an A\ o D n -\ magnonic 
structure. This fact made it natural to generalise the A46J + 012 and A4io,n +012 TBA 
systems in a similar manner, by 'nesting' an A2 o -D n _i magnonic system on top of the 
two ajj soliton nodes. Graphically the result is shown in figure |||. 

Working first at the level of associated sets of functional equations called Y-systems, 
checks on periodicities and central charges were used to fix the details. We finally 
converged to the proposal of appendix |B.2| (Case D), which covers the family of models 
A44 n+ 6,4n+7 + 0i2- This gave a clear signal that the general 0i2 perturbations of unitary 
minimal models share the 'modulo 4' property of the Z^r-related systems found in |5q| , 
an observation which enabled us to obtain the remaining cases by replacing the A2<>D n _i 



part with the other families of systems from [56]. Full TBA equations were then obtained 
by Fourier transforming the Y-systems. 

Once it was seen that the systems in [^] were playing a central role, it was possible, 
following similar reasoning, to conjecture the TBAs for the 021 perturbations. In these 
cases elements of the S-matrix were first used to fix the 'diagonal' part of the TBA 
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Bj B 2 S 
Figure |: Sine-Gordon TBA at ^ = \ 



B l B 3 K 2 



Figure §: + ■ 




Figure |: Sine-Gordon TBA at ^ = f 




Bi B 3 K 2 Kj 

Figure S 7Wio,n + 012 



• • • 

B i B 2 S 



Figure |: Sine-Gordon TBA at = ^4, 




Figure 0: 7W 4 „+6,4n+7 + ^12 



equations, with the systems of ref. [^] then being adapted to describe the magnonic 
partQ. The Y-systems were then derived by using standard methods, paying attention 
to the pole structure of the kernel functions (see appendix A.l and, for example, fl42| ). 
As should be clear from this discussion, the process was one of educated guesswork; 
however the detailed checks that we report below leave us in no doubt that the final 
systems are correct. 

"The two 0(iV)-related set of equations, i.e. those connected to the Zjv TBAs with N even, are very 
similar to the equations for the perturbed 0(N)s/0(N — 1)3 coset models proposed by Paul Fendley 
in | p7[ . The main difference between the two systems can be, naively speaking, traced back to the 
dissimilarity in the neutral bound-state sector. 
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5 TBA conjectures: 2 i perturbations of unitary minimal 
models 



In this section we shall be concerned with thermal perturbations of critical Potts models 
at points where the central charges of the unperturbed models coincide with those of 
the unitary minimal models M PtP +\, p > 3, which means that we set y/q = 2 cos ^pjy > 

and A = |(p— l)/(p+l) , £, = p ■ As explained earlier, the fact that the related minimal 
models are unitary means that we will equally be describing the behaviour of the ground 
states of perturbations of minimal models, the perturbing operator being foi hi all cases. 

For p = 3, 4 and 5, A < 1 and the spectrum consists of the fundamental kink alone. 
TBA systems are already known for these cases - the first is the (trivial) example of 
the thermally-perturbed Ising model, the second was treated in ]3q| , and the third is 
the thermally-perturbed 3-state Potts model, for which the TBA was written in |3S[] . 
Therefore we will suppose that p > 6, which ensures that the bulk spectrum has one 
kink, and one breather M. The form of the 'non-magnonic' part of the TBA is then 
fixed almost completely by the results of §||: 

e B {0) = Rm B cosh6- ^ $ Bc * L c {6) , 

c£{B,K} 

£ k(Q) = RmxcoshO— $Kc * L c {6) + (magnonic terms) , (5-1) 

ce{B,K} 

with the effective central charge given by 

O POO 

c (r) = — / dO m c R cosh dL c (0) , (r = m K R) . (5.2) 

The kernel functions involving diagonal S-matrix elements are: 

W#) = -i^s ^S BB {9) , $kb(0) = <5> BK {e) = -i^ \nS KB {e) . (5.3) 
at) dv 

The last term in ([O]) indicates the presence of the extra contributions from the as- 
yet unknown number of magnons. In addition, equations involving magnons and kink 
pseudoenergies are required. It is at this stage that our conjectures begin, with most of 
the deduction proceeding, as explained in the last section, by analogy with the models 
studied in |3^, [56|. We start with the kernel <&kk(0)- Since the scattering of two kinks 
is not generally diagonal, this is not expected to be a simple logarithmic derivative of 
an S-matrix element. However, if we set 



1 1 

2 ~ 2A 



(5.4) 



then we observe that the mass relation 

tub = 2 cos(7rp)mK (5-5) 
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is accompanied by the following S-matrix identity: 



SaB{0)=SaK(8 + m P )S a K{e-mp) , (a = B) (5.6) 

or, equivalently 

*aB(P) = **K(0 + ™P) + **K(P-i<irp) , (0^O,a = B). (5.7) 

The constraint 9 ^ is due to the fact that &bk{6) has poles at = ±iirp. Care is 
needed when the identity is integrated, either in Fourier transforms or in convolutions. 
Relations between the kernels such as ( [5.5D and (|5.7| ), involving every pair of particle 
species, are needed in order to convert a TBA equation into a Y-system, and so we 
shall assume that ( |5.7| ) holds for a = K too. Assuming also that &kk{6) is free of 
singularities in the strip |Sjm0| < np, we obtain 

W) = x„ • »™w . *<»> = 2pcMh 1 (( , A , ) • (5-8) 

For the reasons outlined in the last section, we looked to the TBA systems which had 
previously arisen in the context of self-dual perturbations of Z^r-symmetric conformal 
field theories for the remaining equations. It turns out that for the <foi perturbations 
of the unitary minimal models A4 PjP +i with p > 5, the equations proposed in |56j can 
be adapted to provide an appropriate set of 'magnonic' equations. The equations in 
[ |56| 1 were obtained via a 'doubling' of the sine-Gordon TBA systems at certain special 
points (3 2 = 32vr/iV, £ = l/(JV/4 - 1). 

The recipe goes as follows. Consider first the continued- fraction expansion of 

(5,) 

For the unitary models this results in four families of cases: 

A) p = An + 2 : 6p ' 

B) p = 4n + 3 : 6p 



1 + 



i 

n-l+i 
1 



1 



n-l+T 



2 

C) p = 4n + 4 : 6p = — — — (5.10) 



D) p = An + 5 : Qp 



1 + 



matching the four families of cases seen in |j5q| . The second 'nested' continued- fraction 
decomposition of 6p in ( |5,10D , i.e. the decomposition of pi = 1/(6/9) — 1, matches the 
special values £ = 1/ (N/4— 1) observed in connection with the Z^r models, if we identify 
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N = p — 1. The fact that the match is at a nested level reflects the idea that these extra 
pseudoenergies are supposed to be of magnonic type, and for this reason we also trade 
the driving terms of |5{| for 

^ a) (0) = 5 al i> lj *L K (9) , (5.11) 

and set 



(magnonic terms) = iplj * Lj (9) , (5-12) 

3 

where j = 1, . . . ,n for cases A, B and C and j = 1, ... ,n + 1 for case D. Finally, 
the variables 9, 9' used in [56] should be rescaled. This can again be interpreted as a 
consequence of the fact that the match is at a nested level and the role of 1/h = 2/{N—A) 
of [56] is now played by the quantity 1/g = pp2 = 4p/ (p— 5) = 2p/h (i.e. g = 3(p— 1)/2). 
This is equivalent to a replacing the kernels <f>^ {9) , tpf^ (9) and 4% (9) defined in 
appendix A of |56| with Pi 



In appendix A the resulting set of TBA equations is written explicitly, while in §^ we 
report some analytical and numerical evidence for their correctness. 

Before that, we will give the explicit forms of the TBA systems for the first four 
cases, p = 6, 7, 8, 9. These are the first unitary models for which the extra particle 
B = Bi is expected to appear, and at the same time their magnonic structures are 
still relatively simple. In each case, the TBA system is conjectured to describe the 
perturbation of A4 PiP+ i by (f>2i ■ 

A) p = 6 (n=l, A = if): 

eb{9) = Rm B cosh9 - ^ ® Bc *L c (9) 

ce{B,K} 

e K (9) = Rm K cosh9- ^ $ Kc * L c {9) - <\> 2 * L w (9) 

ce{B,K} 

6 

e (7)(0) = -^11^^2*^(9)- 5^2* L K (9) ( 7 =1,. ..,6) (5.15) 

0=1 



"Notice that the fact that the kernels involved in different 'nested zones' just differ by a rescaling of 
the arguments of the type: — > 6/(^1^2) — > 6/(^1^2^3) • • ■> with 

* = « ' 6 = SEPTET (5 ' 13) 

is typical of the a^'-related systems |38| (see also fjijl). 
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In this case there are six magnons and zj^ 6 ' is the incidence matrix of the Eq Dynkin 
diagram with nodes labelled as in figure 

B) p = 7 (n=l, A 



9). 



eb(0) 



RrriBCosh.9- ^ 3>Bc*L c (9) 

cG{B,K} 

Rm K cosh 9- ^ § Kc * L c (9) - fa * L m (9) 



c£{B,K} 



( 7 = 0,..., 4) (5.16) 

/3=0 

Here, there are five auxiliary pseudoenergies and l]^ is the incidence matrix of the A§ 
Dynkin diagram with nodes labelled as in figure ^. 

C)p = 8(n=l, X = l): 



eb{9) 


= Rm B cosh 9 - ^ &b c *L c (9) 




ce{B,K} 


ek{9) 


= RmK cosh 9— §Kc * L c {9) — 03 




cG{S,K} 




-03*(L K (0) +j L( 6 )(0)+L( 4 )(0))-0 4 


8^(9) 


= 2 *(if (5) (fl)-L (1) (#)) 




= fo*(K (3 \9)+K^{9) + K^(9)) 


e (6) (0) 


= 2 *id 5) (6<) 


8^(9) 


= 2 *(K (5) (#)-L (2) (0)) 


8^(9) 


= -0 3 * {L^{9) + L^{9)) - 4 * L^{9) 


where K^(9) 


= ln(l + e £(c) W). 



(5.17) 



D) p = 9 (m 

8 B (6 
8 K (6 

eW(6 



1,A = |): 



i?ms cosh 0- ^ §Bc*L c (9) 

c£{B,K} 

Rm K cosh9- ^ * L c (0) + V * + L (3) (#) 

ce{B,K] 

2 

2 ft'V * lW ^ + <W * L K {9) ( 7 = 1, 2) 

/3=1 
4 

2 ft'V * W + <W * (7 = 3, 4) 

0=3 



(5.18) 
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where tp(9) = ipn(0) = ^22 (#)• in this case there are four auxiliary pseudoenergies and 



is the incidence matrix of the A2 Dynkin diagram with nodes labelled as in figure 
It was noted in [f^] that this TBA can be mapped into a particular case of the 



GVrelated Y-systems of [43|. Notice also that the pseudoparticle part factorises into a 
pair of ^2-type TBAs while a single one was found for the model studied in |3f 

This can be explained by observing that the vacuum incidence diagram (k = 8 in 
is (orbifold) equivalent to the product of two "hard hexagon" diagrams (k = 3 in 
Thus also in this case a simple variant of the analysis of |35| should lead to a more 
rigorous derivation of the system, though this remains to be done. 




B K 




Figure @: Mq^ 



Figure |: M 7 ^ + ■ 





B K 



B K 



Figure 0: + 4> 2 i 



Figure pi: M 



9,10 



6 TBA conjectures: A > 3/2 

We were also able to construct the TBA equations describing the ground-state energies 
of the unitary A4 PtP +\ + 4>i2 models. As in the cases related to 4>2i , we split the set of 
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unitary minimal models into four families, which we label A. . . D. For p > 7, the spec- 
trum consists of just four particles (two neutral particles B\ and B% and two kinks K\ 
and K<i). Here the magnonic sector resembles again the Z^r-related structure described 
in A pictorial representation of the TBA systems for the family ^4^+3,4^+4 is 

given in figure [ll]. Note that at n = 0, what for n > had been a magnonic £g-type 
sector in the TBA naturally generates the extra six neutral particles needed to recon- 
struct the full E'g-related mass spectrum of the 0i2-perturbed Ising model. The full sets 



6 




B l B 3 K 2 K l B l B 3 B 5 B 7 B 8 B 6 B 2 



Fi gure 11: Diagrammatic representation for the A / (4n+3,4n+4 + <f>i2 TBA systems 

of TBA equations and Y-systems are given in appendix [B|. 

We shall illustrate these proposals with the simplest new system, corresponding to 
the model M^fi + <j>i2- For this case, A = |g and so, referring back to §[|, we expect 
the spectrum to contain six particles, B±, B 2 , B3, B§, K\ and K 2 . In addition, we 
conjecture that two magnonic pseudoenergies, and e^ 2 \ should be introduced to 
deal with the non-diagonal nature of the kink scattering. The equations are: 

e b {6) = Rm bC oshe-Y J ^bc*L c {9)-{5 btKl +5 bt K 2 )4>*L {1) {0) 

c 

e (i)(0) = -J>*(LW(e) + L Kl (0)+L K2 (e)) (6.1) 

£ ( 2 )(e) = -<p*L (1 \9) 

with b, c G {Bi, B 2 , B%, B$, Ki, K 2 }. The kernels are: 

$w{0) = -i^\nS w {e) , ® kb {Q) = $ bk (9) = -i^\nS kb {6) (6.2) 
where b,b' G {B u B 2 , B 3 , B 5 }, k G {K X ,K 2 } and 

*K 1 K l (e) = <i>** Kl B a {e) , <s>k 2 kM = 0*$k 2 bM, (6-3) 
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§K X K % (0) = {6)=4>* QbsK! {0) , 



A diagrammatic representation is shown in figure |12 



B 1 B 3 B 5 K 2 



O 2 

a i 



21 



cosh(21 i 



K 1 B 2 



(6.4) 



Figure 12: Diagrammatic representation for the Ms,6 + 4>i2 TBA system 



Again the ^2-type pseudoparticle sector reflects the relationship between this theory 
and the hard hexagon model 



7 Analytical and numerical checks 

This section gathers together some evidence, analytical and numerical, in favour of the 
conjectured TBA systems. The value of the UV central charge is obtained through 
standard manipulations of the TBA equations [38, which allow c uv = c(0) to be 
expressed as a sum of Rogers dilogarithm functions. The arguments of the dilogarithms 
are expressed in terms of the set of stationary solutions of the TBA at r = and r = oo. 
We did not recognise these sums among any of the standard relations for the Rogers 
function (see for example |3l|), so we relied on numerical work to check that, up to 
p = 13 and to about 15 digit accuracy, we have the expected result 



c uv = c(0) = 1 



6 



p{p + 1) 



(7.1) 



Next, we discuss the corrections to c(0) which make up the UV expansion of c(r). 
A key idea is to find a set of functional relations satisfied by the exponentiated pseu- 
doenergies, called a Y-system. These generally imply a periodicity property for the 
pseudoenergies under a certain imaginary shift in 9, which in turn can be used to de- 
duce information about the r-dependence of c(r). This phenomenon was first noticed 
by Al.B. Zamolodchikov in ||62|| , for Y-systems related to the TBA equations described 



,65], while a general proof also covering the 



in [40 1 and encoded by a single ADE Dynkin diagram. (The periodicity for the A and 
D cases was subsequently proved in J63|, 
E cases was only given recently, in (6£ 

The Y-systems pertinent to 4>2i perturbations are collected together in appendix A. 
We did not attempt to prove their periodicity properties, but numerically verified the 
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following result: 



Y a (e + mP)=Y a (6), P=_-' (7.2) 



Following the arguments of [62], periodicity implies an expansion for c(r) in powers 
of r^l p and fixes the dimension of an odd perturbing operator to 

a -A-.-I = E±?. (,3, 

This value coincides with the conformal dimension of the operator 02i- 

A numerical study of c(r) for the 021 TBA systems also indicates the presence of a 
single irregular 'anti-bulk' term: 

CO 

c(r) =c-Br 2 + ^ c n r n ( 3p " 3 )/ p . (7.4) 

n=l 

This term is related to the bulk free energy of Q2.14 ) by B = —d£/(jTm^) ; it can be 
calculated analytically from the TBA by adapting the method of |38|, f5J. For 9 — > — oo 
we have 

W) ~ -fcg], f ^ f ^ e e + . . . , (a,b£{B, K}) , (7.5) 

\"W \ m B J 



with 



, (1) . ( , 27T. ,7T. ,7T 7T.\ 

$bb = 4 ( sm (y) + sin (^) + sin (x - 3 ) ) • ( 7 - 6 ) 



Arguing as in j38|, 4C], we therefore have 



vrcD« WW " 2vr sin(f ^) 

This matches the result given in [67]. 

Next, the coefficients c n in the expansion ( |7.4p can be checked against the results of 
conformal perturbation theory (CPT). This predicts an expansion 

6R 00 
c CPT (r)^-—E(r,R)=c+y / B m t m , t = -7tt(R/2tt) 2 ~ 2A , (7.8) 

which should match ( [T.4| ) save for the absence of the anti-bulk term —Br 2 . The CPT 
coefficients B m are given in terms of connected correlation functions on the plane. For 
unitary cases, 

24 I rn—l ffi % 

Bm = m-1 I / (g(M)g(gl.gl)---g(3W-l»5w-l))c H f,- U-A » ( 7 - 9 ) 
JTl. J ^•Zfc.ZfcJ 
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with conformal invariance fixing 



Bo = 12 



7 2 (A) 
7 (2A) 



7(x) 



r(x) 
r(i-x) 



(7.10) 



The high-low temperature duality symmetry of the model is reflected in the fact that 
correlators with an odd number of e operators vanish identically (B m = for odd m). 
Comparing ( f7.4| ) with ( |7.8[) confirms that our effective central charge behaves as that of 
a minimal conformal field theory perturbed by 021 (A21 = A21 = )■ On dimensional 
grounds there must be a relation between r and (the mass of the kink if = ifi 
appearing in (5.1) and (|5.2|)) of the form r = Km^ 2 ^, where k is a dimensionless 



constant whose exact value was calculated by other methods in [67]. Table 2 compares 
numerical results from the TBA with the exact formula. 



Model 


k 2 (TBA) 


k 2 (Exact) 


M 6J + (, 


hi 


0.0259899061737 


0.02598990617395 


M 7 ,s + 4 


21 


0.02506361950681 


0.02506361950686 


M 8 ,9 + 4 


21 


0.0242806346097 


0.02428063460990 


■Mg.io + (, 


hi 


0.02362954340274 


0.02362954340277 



Table 2: TBA data versus the exact results of [p7|1 



Similar checks can be performed for the 4>12 perturbations. The relevant functional 
relations are given in appendix B, and we verified that in these cases the periodicities 
implied by the Y-systems are P = g(p+2) • ^ ne resulting conformal dimensions (p — 
2)/(4p + 4) match those of the 4>i2 operator. 

Finally, we turn to the large- R behaviour, restricting attention to the 4>2i systems for 
brevity. The system (5.1) implies the following asymptotic for the ground-state energy 
E(mK : R): 

E(m K ,R)=--?-c(r)~am K 7rK 1 {r), (7.11) 
bR 

where a can be evaluated in terms of the quantities "Yy in appendix C of |]56f| . Setting 
rj = tt/(j) + 1) and n = Int[(p — 2)/4] we have 

a = |1+ ' 



T (i) 

1 n 




8=1 

( sin(2?7) sin((2n + 2)rj) V^ 2 y-r / sin(2^) sin(ry) 

V sin 2 (n?7) / ~_ ^ V sin((i + 3)rj) sin(ir/) 

4 cos 2 r?- 1 = q- 1 (7.12) 
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for cases A,B and C and 



a 




\ T (i) 

\ n 

( 2 sin(77) cos((n + l)rj) \ t-t / sin(2??) sin (77) 

\ sin(nr/) J y sin((i + 3)77) sin(z7y) 

4cos 2 r/- 1 = g- 1 (7.13) 



for case D. The value a = q — 1 matches the prediction of an argument given in |5f 
extending the instanton ideas of |68|, ^39|, that a should be the Perron- Probenius eigen- 
value of the (in this case S^-symmetric) incidence matrix lj$ encoding the number of 
light kinks which join each pair of vacua. Note that the argument is a little formal 
here, since in general the number of vacua is not an integer. For these unitary cases 
this worry can be averted by shifting to the RSOS approach of [j?]]; either way, the TBA 
systems meet our expectations. 



8 Results from the nonlinear integral equation approach 

The TBA is not the only exact technique on the market for the study of finite-size 
effects. Another class of methods is known as the nonlinear integral equation (NLIE) 
approach. These equations generally depend smoothly on their parameters, and so are 
in some senses better-suited to a description of the g-state Potts models. The downside 
is that the full particle content is encoded in a much more implicit way, and so they 
have less to say on the issues of bootstrap closure discussed in Q . 



8.1 Massive flows 

A non- linear integral equation describing general <pi2, 4>2i an d 0i5 perturbations within 
a unified framework was proposed in [12], similar in spirit to the equations related to 
013 perturbations introduced in [7l| . In contrast to the TBA equations discussed 
earlier, in 0] the set of pseudoenergies { e a } is traded for a single unknown function 
f(9), which solves the following equation: 

f(9) = ina-ir sinh6+ [ ip(6-6')ln(l+e f( - e,) )d6'-[ (p(9-9') hul+e^ 9 ')) d9' , (8.1) 

JCi Jc 2 

where r = m^R with mx the mass of the fundamental kink. Given f(9), the effective 
central charge can be evaluated as 

Ccff (r) = ^ (^j smh9ln{l+e m )d6 - J sinhfl ln(l+ e - /(9) ) d9^j . (8.2) 
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The contours C\ and C2 run from —00 to +00, just below and just above the real #-axis, 
and the kernel (f(0) depends continuously on the parameter A as 



, ,_ r e^sinh(§fc)cosh(f(l-|)/c) dfc_ i d 

As before, yfq = 2sin(^A) , while S(0) is the scalar part of the fundamental kink-kink 
scattering amplitudes as given in (A scale-invariant version of this system had 
previously appeared in the context of lattice models; see f72|.) In fl2]|, we showed 



how to use (8.1) to describe perturbations of unitary and non-unitary minimal models, 
tuning the twist parameter a so as to account for the fact that the ground state may be 
generated by a negative-dimension field. For current purposes, we are more interested 
the q-st&te Potts models and their tricritical variants, for which, as discussed earlier, 
the ground state always comes from the identity operator. For these theories, a should 
be related to A as 

« = y-l. (8.4) 



The reader should note that the NLIE ( p.2|) is simpler and, not being restricted to 
special points, of wider applicability than the TBA equations proposed in §[| and §^ 
above. It is also much more appropriate for discussions of the q — > 1 limit relevant 
to percolation. However, it does not encode, in any direct way, information about the 
bound-state content of the associated field theory. 

In Table 3 we compare results from the NLIE and TBA approaches. The agreement 
is very good, ranging from 10 to 14 significant digits. 



Model r 


TBA NLIE 


.Mio,ii + 02i 0.2 
0.3 


0.9328551607515 0.9328551607514 
0.9185246910885 0.9185246910884 


•Mll,12 + 021 0.2 
0.3 


0.9421554960427 0.9421554960424 
0.9279952931412 0.9279952931410 


■Mi2,13 + 021 0.2 
0.3 


0.949321521251 0.949321521253 
0.935305379646 0.935305379648 


A4l 3 ,14 + 021 0.2 

0.3 


0.954961250612733 0.954961250612730 
0.941068318324234 0.941068318324237 


A45,6 + 012 0.2 

0.3 


0.79241547642 0.79241547640 
0.78313881742 0.78313881742 



Table 3: Comparison between TBA and the NLIE proposed in |T!| 
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Temperature 



Potts CFT 



Disordered phase 




Ordered phase 




Tricritical Potts CFT 



Fugacity 



Figure 13: The phase diagram 



8.2 Massless flows from tricritical to critical models 

In addition to the perturbations discussed so far, the concentration of vacancies in a 
tricritical Potts model can be adjusted so as to drive the theory either onto a line of 
first-order transitions, or else down via a massless flow to the corresponding critical 
model. The relevant operator turns out to be cftis, which is always in the spectrum of 
the tricritical models. The full picture is illustrated in figure [l^; it is most convenient to 
use the parameter £ for these flows, related to q by yfq = 2 cos 7r/(£+l) with £ £ [1, oo]. 
The massless flow is then 

^^(0 = 1- (€ + i) 6 (e + 2 ) Ccrit(0 = 1 " WPT) ■ (8 ' 5) 



For £ = p G Z, this is the well-known massless 4>i3 flow from -M p +i jP +2 to M PtP+ \ [ 73|, [74[] 



for which equations of TBA type were proposed in [75]. To describe the flows at general 
values of £, we instead adapt an equation proposed by Zamolodchikov in 111] for the 



study of flows in the imaginary-coupled sine-Gordon model from c = 1 to c = 1. In [37], 
it was shown that the twist parameters in Zamolodchikov's equation could be chosen so 
as to describe flows between both unitary and non-unitary minimal models, revealing a 
striking non-monotonicity of c(r) at intermediate scales in the non-unitary cases. Here, 
we shall propose a further modification to capture the general flows between tricritical 
and critical Potts models. As in [13], introduce two analytic functions fn(0) and 
and couple them together via 

fR(P) = -i r -e e + ma' 

+ ( 0(6 -e')ln(l + e fR(e ">) dO'- [ 0(9 - 6') ln(l + e~ fRie,) ) d9' 
JCi Jc 2 
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+ / X(P- 0') ln (! + e- /L(9 ' } ) dB' - [ X (0- 0') ln(l + e /i(e,) ) dB' 
JCi Jc 2 



h(0) 



T 



nra 



[ <j){e - e')\n{i + e fL{e>) ) de' - [ </>(0-0')M 1 + e ~ /£(0 ' ) ) d0 ' 

JC 2 JCi 



where 



+ f x{0-0')Ml + e~ Me,) ) d9' - [ x{0-0')H^ + e M8/) ) d6' (8.6) 
JC 2 JCi 

"°° e ike smh(k(£- l)f) dk 



x{0) 



2cosh(p)sinh(A;f £) 2vr 



(8.7) 



e ik6 ginh (Z^ 



dk 



2cosh(^)sinh(^) 2vr 



(In yj|] , the kernels were given in terms of a parameter p = £+1 .) We have = 
2 cos 7r/(£+l) with £ G [1, oo], and, for the interpolating Potts flows, a' must be chosen 
as 



a 



(8.9) 



In terms of the solutions to these equations, the effective central charge is 



Ceff(r) 



3ir 
2^2 



f e e hx(l+ef"W) d9 - [ e e \n{l+e- fR{e) )d6 
JCi Jc 2 

+ f e- e ln(l+e /L(e) ) dO - [ e~ e ln(l+e- /L W) dB 
JCi JCi 



UO) 



Figure 14 shows our numerical results. As long as the IR destination has a positive 



central charge, the flows are monotonic, even at points corresponding to non-unitary 
minimal models. This is in sharp contrast to the behaviour of the (^13-perturbed minimal 
models themselves, for which a number of non-monotonic flows of c(r) were exhibited in 
To some extent, this reflects the already-mentioned fact that the ground states of 
Potts and minimal models do not coincide away from unitary points, but it is nonetheless 
surprising that the switch to the identity vacuum state manages to eliminate all of the 
non-monotonicity while cjr remains positive. 

The other notable feature of figure 14 is the behaviour of the flow from c = 
to c = —2, which exhibits a cusp, suggestive of a level-crossing, at an intermediate 



length-scale. In fact, this curve deserves a second glance for another reason: in |7q| , 
Fendley, Saleur and Zamolodchikov pointed out that the effective central charge along 
this particular flow should be protected by supersymmetry, at least for small enough 
values of r, and hence should be identically zero. They also, by an indirect argument, 
predicted that this picture would be changed by a level-crossing at r w 2.95708396. The 
apparent contradiction of their first claim with our results has a neat resolution: one 
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-4 



-3 



-2 



-1 







Figure 14: Flows of c(r) as obtained from the NLIE flOP, plotted against log(r/2). 
Highlighted are the flows from c = 1/2 to c = and from c = to c = - 2. 



has to remember that, just as for the TBA equations discussed earlier, the quantity 
c(r) produced by the NLIE includes an 'anti-bulk' piece, non-analytic in the coupling 
constant, which ensures that at large values of r, c(r) tends to a constant. Thus, just 
as in equations ( k-J ) and (7 



above, we expect 
c(r) = c CPT (r) - i3r 2 



ill 



where c CPT (r) is the physical, unsubtracted quantity, while c(r) is the quantity which 
is found directly from the NLIE. The 'Bulk' constant B was found exactly in [77| ] by 
considering a related theory in a magnetic field; it can also be obtained directly from 
equations ( p.6[ ) and ( |3.10| ), using an argument described for a similar equation in |56|] . 
Either way, the answer is 

3 1 

^ = -27cos(f(e+i)) • (8 - 12) 

Specialising to £ = 1, we see that before the bottom curve of figure [14] is compared with 
the proposals of JF"6|], the quantity ^r 2 should be added to it. Numerically, it is hard to 
do this directly due to instabilities in the equations near £ = 1, so in figure [D] we show 
a sequence of plots of the appropriately- adjusted functions Ccpt 

(r) = c(r) +B(C)r 2 , for 

£ = 1.03 down to £ = 1.005, together with an extrapolated curve for £ = 1. (For the 
extrapolation, numerical data down to £ = 1.001 was used.) Not only does the revised 
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curve meet general super symmetric expectations; the point at which supersymmetry 



was predicted in [76] to be spontaneously broken via a level-crossing is also reproduced. 
We suspect that there is more to be said on this matter, and the nice agreement between 



our results and those of |7g] certainly deserves to be understood at a deeper and more 
analytical level. For now, we leave this for future investigations. 




Fi gure 15: Flows of c^pt 

(r) for £ = 1.03, 1.025, . . . 1.005 as obtained from the NLIE, 
plotted against log(r/2). Also shown is the extrapolated curve for £ = 1 (dotted), and 



the location of the level-crossing predicted in J76J. 



9 Conclusions 

In the past decade there has been a common belief that A n -, D n — and E n - related 
integrable systems should exhibit a certain uniformity of structure. As we remarked in 
j|| , the bulk, S-matrix, description of the real-coupling simply-laced affine Toda field 
theories (and their minimal variants) illustrate this philosophy well: all the scattering 
data can be encoded using the root systems of the associated Lie algebras [Q. A first 
sign of a breakdown in this pattern in other contexts came with the attempt to make 
the coupling imaginary in these same Toda models: the A\ Toda Hamiltonian remains 
Hermitian, turning into the sine-Gordon Hamiltonian, while all the others become com- 
plex. Despite this crucial difference, it was initially thought that the relatively simple 
j4i-related sine-Gordon spectrum might be the first instance of some unified description 
valid for all Lie algebras |79, pQ] . Smirnov's work M had already suggested that matters 
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were likely to be much more complicated, and subsequent work, in particular |yj and 
has confirmed this expectation, to the extent that there is still no complete un- 
derstanding of the bulk scattering theory associated with any imaginary-coupled Toda 
theory beyond the A\ case. 

On side of lattice models and the Bethe ansatz, a longstanding hope has been to 



find generalisations of the construction of Takahashi and Suzuki |58|, who many years 
ago found the correct 'string hypothesis' describing the Bethe ansatz solutions relevant 
to the finite volume corrections-to-scaling of the XX Z model at arbitrary values of the 
anisotropy. The form taken by the Bethe ansatz equations for this model are the Ai 
cases of a general set of ADE-related Bethe ansatz equations, given just in terms of 



Lie algebra parameters in [82]. The symmetry of equations in nature is often broken 
by their solutions, and for three decades all attempts to generalise the construction by 
Takahashi and Suzuki to the other lattice models at arbitrary values of the coupling 
constant have failed. (The difficulties even for the SU(3) case were highlighted in |]8lfl .) 
Building on a relationship with the q-state Potts models, in this paper we have 

(2) 

conjectured and checked TBA equations related to a 2 at points related to the unitary 
series of minimal models perturbed by their <p2i and 4>i2 operators. Although we have 
proposed a continued-fraction decomposition, a la Takahashi-Suzuki, that governs these 
systems, we are currently unable to extend the analysis to the other rational points. Yet, 
the models described, being unitary, are physically the most interesting, and we hope 
that our work will motivate a more rigorous derivation of these equations. In this 



respect we feel that the functional approach described in, for example, [83, 84 1, and 
already applied by J. Suzuki to particular examples of $12 perturbations in [|l7], 18], is 
likely to be the most effective. For the Potts models, the TBA equations that we have 
been able to find serve to confirm the mass spectrum found by bootstrap techniques in 
at least at the points = 2 cos(-7r/(p+l)). 

Finally, we mention two further open problems. First, it would be nice, along the 
lines of |8^, |86|, |87], 61], to adapt both the NLIE and the TBA equations to describe 
excited states. (The recent paper [^9| discusses TBA-like equations for excited states in 
one particular (^21-related model.) In particular, we note that the ground-state NLIE 
of § |8.1| works without problems even in the region A € (|,3] where we encountered 
problems closing the bootstrap in Q. Having control over the full finite-size spectrum 
using the NLIE technique should help us to resolve this open question. Second, we have 
introduced many new sets of TBA equations and Y-systems in this paper. It seems 
important to elucidate the associated character identities^ to find the T-systems |84j] 
and to give to the T functions a spectral interpretation in the spirit of the ODE/IM 



correspondence [12|. Relations with (generalised) quantum KdV equations pOL B3, 91, 



|93j] and with perturbed boundary conformal field theory |90|, |9lJ might also be explored. 
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A TBA and Y-systems (^21 perturbations) 

In this appendix we list the TBA equations and Y-systems for the Ai PlP +x + 4>2i models. 
As mentioned in §|| above, the continued-fraction expansion of the parameter 6p = 
J allows us to distinguish four families of systems: 



A) p = An + 2 : 6p 

B) p = 4n + 3 : 6p 



1 



1 + 



n-l+i 
1 



1 + 



C) p = 4n + A : 6p = — , (A.l) 



ra-H — K 



D) p = 4n + 5 : 6p 



1 + i 

n 



A.l The kernels 

The kernels functions related to the diagonal S-matrix elements were defined in (|5.3| ), 
and we do not repeat them here. 

There are also kernels associated with the magnonic pseudoenergies, which were 
given in equation ( 5.14| ) of the main text in terms of those of [56]. For completeness, we 



give them here explicitly. By analogy with a notation used for the affine Toda theories, 
define the blocks 



Then 

^= II UH 1 -^' (j,k = l...n-l) , (A.3) 



j+k-i 



step 2 
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and 

h/2+k-l 

S kn = Sfc,n+i = (-1)" II U}> (fc = l---n-l) > (A.4) 

h/2 — k+l 
step 2 

where h = {p — 5)/2. Some of the magnonic kernels can now be expressed in terms of 
the logarithmic derivatives of the S functions: 

<j> jk {9) = -i-^\TxS jk {9/2 P ) , i> jk (e) = -i-^lnf jk (e/2p). (A.5) 

In the definition of ipj k , the function Tj^ is obtained by replacing each block {x} in 
( |A.3j )- (|A.4|) by (x). In ( A.5 ), j and /c run from 1 to n—1, with n depending on p as in 

For the D-type models h = 2n and the definition can be extended immediately to 
cover the remaining cases when both j and k take the values n or n+1. The remaining 
functions needed are given by 

2n-3 2n — 1 

neven : S nn = S n +i, n +l = ] [ { | } , S n<n+ i = ]^[ {^} , 

1=1 i=3 
step 4 step 4 

2n— 1 2n-3 

n odd : S nn = § n +l,n+l = ]~][ {^} j Sn,n+1 = — {^} • 

(=1 1=3 
step 4 step 4 

Otherwise, the associated kernels are more elaborate. Define an integer t by p = 
4n+t + 1, and then set g = 3(p — 1) /2 and 

= (A ' 6) 

This function has the property that 
Then the kernels in the TBA are given by 

<Pnn(0) = Xt*<f>n,n-l{9) > 4>nn(0) = Xt*lpn,n-l(0) ■ (A.7) 

(For p<6, 4> n ,n-i and ip n>n -i are not defined and we set 0nn = VVm = O.) Finally, depending 
on the continued- fraction expansion of 2/g, a number of extra kernels are needed. These 
are 

9 



cosh g0 

2g 8g cosh ^(9 

^ = c^ih2^ ^ = 3(4cosh 2 ^-3) 

h , , 8g cosh 
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This completes the definition of the kernels as needed for the TBA equations. How- 
ever, when deriving Y-systems from the TBA (as we did for these 021 cases) it is also 
important to know the precise locations of the poles in the kernel functions. This is 
because the derivation involves complex shifts in the rapidity 9, and care must be taken 
when poles cross the contours of convolution integrals, as extra terms are generated 
which enter into the Y-systems in a crucial way. 

In most cases these pole locations are easily read off from the explicit formulae, but 
so far we have only given the kernel &kk m terms of an integral representation ( |5.7| ). 
Here we show how an alternative product formula can be obtained. 

Let us define the function B(9) such that 

B(9 + ip)B(0-ip)=S KB (9) (A.9) 

and set 

$kk = -i-^lnB(9) . (A.10) 
Inverting QA.9| ) using Fourier transform we find 

where fKB(k) is the Fourier transform of 

f KB {9)=\nS KB {9). (A.12) 
In order to get a convergent product representation for B{9) let us set 

± _ T(l jf »(fl/7r T ia)/2) T(l T Hd/v jf i(l - a))/2) 

> T{Ti(9/7r±ia)/2) T{Ti{9/-K±i{l-a))/2) 1 ' 

and write 

a+(9,a 1 )a+(9,a 2 ) 
a (9,a 1 )a (9,0,2) 

withai = 1/2+1/2A and a 2 = 1/6+1/2A. Note that cr + (9, a) ( <r~(0, a)) has only a finite 
number of zeroes and poles in the upper (lower) half plane Ssm(9) > (Qm(9) < 0). 
Expanding the cosh function in the denominator in (A.ll) we can formally write 

n=0 ^ ' 

°° / roo \ 

x ]Jexp (-l)( n+1) / dke^ e - i ^ {2n+1 ^ k {a-{k,a 1 )+d-{k,a 2 ))) 
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and hence 



By) = n 

n=0 

oo 

- n 



n=0 



a + (9 + ifm(An + 1), a{) a + (9 + ipTi{An + 1), a 2 ) 
a+{9 + ip-n{An + 3), ai) a+{9 + ip7r(4n + 3), a 2 ) 

<7~(# — ipir(An + 3), ai) cr~(# — i/)7r(4n + 3), 02) 
<7~(0 — ipir(An + 1), ai) cr~(6> — ip7r(4n + 1), 02) 



To derive the Y-systems, the following mass relations were also important: 



7T p— 5 

6 p-1 



Illfi = 2 COS 

2cos (f iT~r) m £ = + 2 cos (f^rr) ™x 



A. 2 General notation for Y-systems 

We define 

Y{9) = e 

and introduce the shorthand notations 



e(0) 



A(±m,. 



m 

, ±n m ) = JI (1 + Y(9 + iir%)) (1 + Y{6 - iw%)) 



and 



C(±ni, . . . ,±n r , 



n 



i=i 



Y(e + iir%] 



1 + 



Y(9-iir%) 



(A.15) 



(A.16) 



(A.17) 

(A.18) 
(A.19) 



(A.20) 



with rij > 0. Just a single factor appears on the RHS for entries with the ± omitted, so 
that, for example, 



A(0, ±n 2 ) = (1 + Y{9)) (1 + Y(9 + ivrf )) (l + Y(9 - 



(A.21) 



In practice the functions e, Y, y and so on appear with indices to show which pseu- 
doenergy is involved. In all of the above definitions, these indices take the same values 
in all factors. 

We shall also set H = 6|£— 1| where, as before, c = 1 — 6/£(£+l) with £ > on the 
critical branch, and £ < on the tricritical branch. For the unitary minimal models 
Mr, p+ i this translates as H = 6(p— 1) for the ^21 perturbations, and H = Q(p+2) for 
the 0i 2 perturbations. 
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A. 3 TBA equations and Y-systems for M. PtP +i + 2 i 

In the figures below, we supplement the four families of $21 -related TBA equations and 
Y-systems with four sets of graphical representations. These graphs give a rough idea 
of the structure of the TBA equations, and also fix the labelling conventions. On the 
magnonic nodes, to any pair of numbers (i, a) corresponds, in the TBA equations, quan- 
tities labelled with a lower index i and an upper index (a): ef , lJ q) = ln(l + e~ £ i ). 



Case A, (p=4n+2 , n > 2): 




eb(0) = Rm B cosh (9 - ^ & Bc *L c (9) 

ce{B,K} 

n 

e K {6) = Rm K coshe- ^ <f> Kc * L e {0) + £ i> v * Lf (0) 

ce{B,K} 3=1 
n 

efXe) = 5 al ^ u *L K {9)-Y J [&3*Lf\e)-^ i3 *Lf\e) 

3=1 

6 

Y. l[ S ] fo* L n\ e ) (i = l,...,n; «=1,2; a = 3-a) 
0=1 

6 

4 7) W = -E^^^fW (7 = 3,4,5,6) 

/9=1 

The Y-system is 

3>b(±(p + 3)) = A x (0,±8)Af ) (±4)Al 2) (0) 

n-l 

y K (±(p-5)) = A B (0)£W(0)l\4\±(4n-4i-3)) 
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n-1 jlAn-l] 

yf\±4) = [£x(0)]^ 1<5Ql 4 a) (0)n(Ai a) (0)) fcJ 

k=l 

A® (0) A<f> (0) A(f) (±1) A(f) (±2)A« (±3) 

A(f) (O)A(f ) (O)A(f ) (±1) A( 4 ) (±2) A( 2 ) (±3) 
(with j = 1, . . . , n — 1; a = 1, 2 and a = 3 — a 



<5, , ,n-l<5 Q! 2 



yW(±i) = a^(o) 



,(7) 



n [^(o) 



l /3 7 



( 7 = 1,2,..., 6) (A.22) 



0=1 



Case B, (p=4n+3, n > 2): 



O (n.4) 

(1,2) (2,2) (n-1,2) 

O — O 0=£O ( n -2) 



(i,D 



6 (n.0) 



O toO (n,D 



(2,1) (n-1,1) 



O (".3) 



B K 



e B (0) 
e K (0) 



RrriB cosh - V] & B c*L c (9) 

c£{B,K} 



RrriK cosh 0- *if c *L c (0) + ^Vii*4 W 

ce{B,i^} j=i 

n 

^i^ii * [<^*4 a) (0) - i>ij*LfXe) 

3=1 
4 

-^nJ^^Vi^W (i = l,...,n; a=l,2; a = 3- 

13=0 

4 

-Eft 5l ^* L ^W (7 = 0,3,4) 



■a) 



(A.23) 
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The Y-system is 

y B (±(p + 3)) = A K (0,±8)A«(±4)A? 3 (0) 



n-l 



y K (±(p-5)) = A B (0)£« (0) n cf\±{An-Ai- 2)) 



n-l 



A a )fn\Y kj 



fe=l 



A(°)(0)A( 3 )(0)A( 1 )(±2; 



^7,n — l<5a2 



A(°)(0)A( 4 )(0)A( 2 )(±2)_ 
(with j = 1, . . . , n — 1; a = 1, 2 and a = 3 — a) 

^ } (±2) = [A^CO)] 71 72 [J [4«(0)] " ,(7 = 0,. ..,4) (A.24) 

/3=0 



Case C, (p=4n+4 , n > 2): 



(2,1) (2,2) (2,n-l) 

o — o- 




B K 



e B (0) 
e K (0) 



Rttib cosh 9 — V] <£ Bc *L c (6>) 

n 

.Rm^cosh^- ^ $ Xc *L c (6') + ^Vii*4 1) ^) 
ce{B,A'} j=i 



n 

^ * l k {6) -J2 [^*Lf\e) - ^*Lf\e) 



$in(>a2 



<M4 4) W + 4 6) («0) + 04*4 3) (^) + <fc*4 5) W 
<M4 3) (0) + 4 6) (*)) + &*4 4) (*) + &*4 5) (0) 
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<t> 2 *(K^(9) + K^(0) + K^(9)) 
fo*K^(9) 

J> 2 *(KW(9)-LW(9)) (A.25) 
with % = l,...,n, a=l,2, a = 3-a, and K^ c \9) = ln(l + e e<c) W) as in (|T7]) above. 



4 3) W 

e^(9) 

4 6) W 

e^{9) 



The Y-system is 

y B (±( P + 3)) 
y K (±( P -5)) 

3>j a) (±4) 



y^(±3) 

^ 2) (±3) 

3^„ 4) (±l) 
3^(±l) 

yW(±i) 



A i ,(0,±8)AS 1) (±4)A[ 2) (0) 



n-1 



A J? (0)4 1 )(0)[]/:f ) (±(4r l -4i-l)) 



i=i 



n-1 



[£K(0)]^ l5al 4 a) (0) JJ (A< a) (0) 
Af(0)AW(±l) 



kj 



k=l 



An 4) (0)An 2) (±l)_ 

(with j = 1, . . . , n — 1; a = 1, 2 and a = 3 — a) 



A 



(2 



1 (o)4 4 )(o)4 6 )(o)4 5 )(±i)4 3 )(±2) 
1 (o)4 3) (o)4 6) (o)4 5) (±i)4 4) (±2) 

(O)Af(O) 
4 2) (0)A( 5 )(0) 
A(f)(0)A( 6 )(0)A( 4 )(0) 
A(f)(0) 



(A.26) 



Case D, (p=4n+5, n > 2): 

eb{9) = Rm B cosh9 - ^ $ Bc * L c (9) 

ce{B,K} 

n+1 

e K {9) = Rm K cosh9- <f> Kc * L c (9) + £ ^ y * Lf (9) 

ce{B,K} 3=1 
n+1 

et\e) = 5 al ^ u *L K {9)-Y\^Lf\9)-^ i3 *Lf\9)\ (A.27) 

3=1 

with i = 1, . . . , n + 1, a=l, 2, and a = 3— a. 
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(n+1,2) 




B K 



The Y-system is 

y B (±(p + 3)) = A i ,(0,±8)AS 1) (±4)AS 2) (0) 

re-l 



n+l _ _ jPn+l] 



^ (Q) (±4) = [£ K {0)] s ^£f\0)H [Af } (0) 



i=i 

(with j = 1, . . . , n + 1 , a = 1, 2 and a = 3 — a) (A. 28) 
A. 4 Exceptional 2 i Y-sy stems 

The well-known |3^] Y-systems for ^3,4 + (£21 and + ^21 are 

3M±6) = 1 (A.29) 

and 

2 

^■(±8) = J] [A,(0)]^ 2] (j = 1,2) (A.30) 

k=l 

respectively, with 3^1 = 3^2 = 3^ in the second system. 

For «M4,5 + </>2i, a Y-system can be derived from the TBA of [j35f|. In the current 
notation, it can be written as 

y (2) (±l) = £ (1) (0) 
yW(±l) = CW(0)/C K (0) 
■v f+f(\ ^(0) £^(±4)^(^1,^5) 

(The periodicity implied by this system is P = 32/ H = 16/9, and the resulting confor- 
mal dimension A = l — 1/P = 7/16 matches that of the perturbing operator </> 2 i-) 
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The remaining exceptional 021 Y-systems can be derived from the TBA equations 
given in §[| above, and are as follows: 
M&j + 02i : 

y B (±9) = A x (0,±8)A (1) (±1,±7)A (3) (0,±2,±6)A( 5 )(±1,±3,±5) 

x A( 6 )(0,±4)A( 4 )(±2,±4)A( 2 )(±3) 
y K (±l) = A B (0)£«(0) 



6 ; 

yW(±l) = [£k(0)} 5j1 II ( £(/3) (0)) 

0=1 



[Be] 

01 



(7 = 1,..., 6) 



(A.32) 



-^7,8 + 021 : 

y B (±10) = A^(0,±8)A( 4 )(0)A( 2 )(±2)AW(0,±4)AW(±2,±6)A( 3 )(±4) 
y K (±2) = A B (0)£«(0) 



l l3l 



(3=0 



(7 = 0,..., 4) 



M 8 .9 + 



8,9 T 921 



y B (±n) = A^(0,±8)A«(±3,±5)A( 2 )(±l)A( 3 )(±4)A( 4 )(0) 

y K (±3) = A B (0)£«(0) 

>(±3) = £^(0)£( 6 )(0)£( 4 )(0)£( 5 )(±1)£( 3 )(±2) 

>(±3) = £ (6) (0)£ (3) (0)£ (5) (±1)£ (4) (±2) 

)(±1) = A( 5 )(0)£ (1) (0) 

>(±1) = A( 5 )(0)£ (2) (0) 

)(±i) = a( 6 '(o)aW(o)a^(o) 

)(±1) = A®(0) 



y {1 

y(2 

y(6 



A*9.10 + 



9,10 T 921 • 



y B (±12) = A^(0,±8)AW(0)A^(0)A( 1 )(±4)A( 3 )(±4) 

yjc(±4) = A B (0)£«(0)£( 3 )(0) 

y«(±4) = £^(0)£( 2 )(0) 

^(±4) = £ K (0)£^(0) 

y (2) (±4) = £«(0) 

y (4) (±4) = £®(0) 



(A.33) 



(A.34) 



(A.35) 



B TBA and Y-systems (0i2 perturbations) 

The mass spectrum for the theory M. PjP+ i + 0i2 for p = 3, 4, 5 and 6 contains 8, 7, 6 and 
6 particle types respectively. In the region p > 10 there are instead only four particles: 
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two breathers ( B\ and B 3 ) and two kinks ( K\ and K 2 ), with masses 

TdBx = 2m COS (£i7T + 4^27r) 

rag, = 4m cos (^i7r + 4^27r) cos (^i7r) (B.l) 

= 2m cos (^tt) 
m^ = m 

where & = (p- 2)/(12 + 6p) = (p - 2)/# and £ 2 = 2/(6 + 3p) = 4/fl". The four masses 
also satisfy the following 'fusion' relation 

2 cos (£ 2 7r) m Kl = m K2 

2 cos (£ 2 7r) = m Kl +m B3 (B.2) 

2 cos (£i7r) m fi3 = m Bl + 2 cos (£i?r - £2^) wiif 2 

2 cos (£i7r) rriflj = mg 3 

We shall introduce a tower of magnonic pseudoenergies determined by the ratio 

p = 6/6 = (p-2)/4. (B.3) 

As for the <foi perturbations, there are four distinct families of models, determined by 
a continued fraction decomposition: 

> + ; 
>4 



A) 


P 


= 4n + 3 


P = 


(n 


5) 


P 


= 4n + 4 


P = 


(n 


5) 


P 


= 4n + 5 


P = 


(n 


D) 


P 


= 4n + 6 


P = 


n 



1+1 



B.l The kernels 

First we need the kernels involving the breathers and kinks. These are simply defined 



as 



<S>ab(0) = <S>ba(0) = -i-^\nS ab (9) , (a = B 1 ,B 3 ,K 1 ,K 2 ; b = B u B 3 ) (B.4) 



where the quantities S a b(9) were defined in §3.3. The kink-kink kernels are defined, as 



for the <?!>2i -related cases in §|5], to match the mass fusion relation. 

^K 2 K 1 = §K^K 2 (0) = W* §B z Ki (0) ~ w(9) 
®KiKi(9) = z*<5> K2Kl {6) 

$k 2 kM = $K lKl (e) + <!> B3Kl (9) (B.5) 
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with 

H , n s H cosh(H6/24) _ , 

Z{e) = 8cosHH9/8) ' W ^ = iVl cosh(^/8) " (R6) 
Some of the extra kernels needed are defined in term of S-matrix elements of the A2, 
A§ and Eq 'purely-elastic' scattering theories. In general these S-matrix elements can 
be written as (78|, ||] 

s\f(o)= n to. ( b - ? ) 



and from these we can define the functions 



where 5 is A2, A§ or an d A^J is a set of rational numbers depending on i, j and g, 
and the blocks (x) and {x} were defined in ( |A.2| ). Then the kernels needed in the TBA 
are defined as 

^(0) = -ip^((p + 2)l) , ^(9) = -i±Tj?((p + 2)l), (B.9) 

and 

4f* ] W = 4f ] (° ~ ™I H ) + ^ (° + i7r /#) • (B.10) 



B.2 TBA equations and Y-systems for M. p . p+ \ + 



>V2 



As in appendix A. 3, we illustrate our proposals with a set of graphs, which give a rough 
idea of the structure of the TBA equations, and also fix the labelling convention. Where 
convenient, we also refer to kink-related quantities using a single label K a , rather than 
the pair (0, a). For the Y-systems we use the notation defined in § A. 2 with H = 6(p+2). 

Case A, (p=4n+3, n > 1): 

e Ba {0) = Rm Ba cosh 9- ®B a , c *L c (6) (a = 1,2) 

ce{B u B 2 ,Ki,K 2 } 

ef\0) = 5 i0 (i2mK a cosh(9- £ $ Ka4 * L d (9) + £ ^ ] * Lf (9) 

\ de{B 1 ,B 2 ,K 1 ,K 2 } /3=1 

2 / n-1 \ 

,8=1 V fc=0 / 

6 

+ S j>n ^J2^aS 6] * L ^H0) (j = 0,...,n-l,a = l,2) 
5=1 

4 7) w = -E^ 6l *4 5 He)+^ 6l *4 1 liW+#*4 2 2iW (B.11) 



5=1 
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(11,6) 



Bj B 3 K 2 

(0,2) 



(n-1,1) 




Ki 

(0,1) 



with 7 = 1, . . . , 6. 
The Y-system is 

y Bl (±( P -2)) 
y B3 (±(p-2)) 



jf } (±4) 



x 



X 



yW(±i) = 



A B3 (0) 

n-l 

A Bl (0)A( 2 )(0)nAf ) (±(4n-4.-3)) 
i=0 

n-l ,[A n ] 

[A B3 (0)]^^ A f (0)n(4 Q) (0)) " 

fc=0 

4 4) (o)4 6) (o)4 5) (±i)4 3) (±2)4 } (±3)] 5j - i5q1 
4 3 ) (o)4 6) (o)4 5) (±i)4 4) (±2)4 2) (±s/ 

(with j = 0, . . . ,n — 1; a = 1, 2 and a = 3 — a 

1 ;l B 6l 



fij,n — 1<5 Q 2 



4 7 ii(o) 



5 7 i+<5 7 2 

0= 



6 Z [B 6 1 

n[ A ^(°)]" 7 (7 = 1,... 6) (B.12) 

13=1 



Case B, (p=4n+4, n > 1): 

e Ba {0) = Rm Ba cosh 9- $B a ,c*L c (9) (a = 1,2) 

ce{Bi,B2,Ki,Jf 2 } 



= 5 i0 [#m^ Q cosh0- 

\ de{Bi,B 2 ,-f!:i 
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(n,4) (n,2) (n,0) (n,l) (n,3) 

O O Q O 



(n-l,2)( 



O (n-l,D 



(2,2) Q 



(1,2) 



_! 



O (2,1) 

O (Li) 



B j B 3 K 2 

(0,2) 



Ki 

(0,1) 



n-l 



p=i 



k=0 



4 7) W = 

with 7 = 0, . . . , 4. 



+ d^Y^Wa * L( nH0) (i = 0,...,n-l;a = l,2) 
5=0 

= - E * 4 5) («o + ^ * l&w + <t 1 * 4 2 4w 

5=0 



The Y-system is 

y Bl (±(p-2)) = A B3 (0) 

y Bs (±(p-2)) = A Bl (0)A( 2 )(0)nAS 2) (±(4n-4i-2)) 



n-l 



i=0 



n-l 



,[A n ] 



^ a) (±4) = [A B3 (0)]^ 5 -Af ) (0)n(4 Q) (0) Vj 

fc=0 

I <5j,n-l<5al 



(B.13) 



4°)(o)4 3 )(o)4 1) (±2)j 
4°)(o)4 4 ) ( o)4 2) (±2)]^- i5q2 

(with j = 0, . . . , n — 1; a = 1, 2 and a = 3 — a) 

c i r 4 ,[A 5 ] 

3^ (±2) = [4^(0)] 71 " II [4^(0)]^ (7 = 0,... 4) (B.14) 

13=0 
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Case C, (p=4n+5 , n > 1): 



(n,6) 



(n-1,2) 




B 1 B 3 K 2 

(0,2) 



ebM = Rm Ba cosh 9- $B a ,c*L c (9) (a = 1,2) 

ce{B u B2,K 1: K2} 

ef\0) = 5 j0 \ Rm Ka cosh 6- £ * Ka4 * L d (9) + * % 



- n 



4 3) (* 



4^ 



d£{B 1 ,B 2 ,K 1 ,K 2 } 
n-l 



0=1 



k=0 



/3=1 

-0 2 *(4 6 )(#) + 4 4 ^) + 
2 *(i4 1} w-4 5) (0)) 

02 * ( 



(j = 0,...,n-l;a = 1,2) 



03 * (lf< 6 >(0) + 4 3 )(0) - L^*)) + 05 * #< 5) (*) + 04 * *< 4) (*) 

03 * (k®(P) + 4 4 )(0) - 4^(0)) + 05 * ^ W + 04 * «^(fl) (B.15) 



where K^ c \6) = ln(l + e e(c> ^^), as in ( 5.17 ) earlier, and the kernels 0j are defined in 
QOP with g = 3(p + 2)/2. 
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The Y-system is 

y Bl (±( P -2)) 
y B3 (±( P -2)) 



X (a) (±4) 



^ n 

V(2 

n 



V(4 
■J 1 n 

y(5 
•s n 

y(6 



(±3) 
(±3) 



(±1) 
(±1) 
(±1) 



fc=0 



'kj 



A B3 (0) 

n— 1 

A Bl (0)A( 2 )(0)nAf ) (±(4n-4i-l)) 

i=0 

ii—i 

[A B3 (0)]^Af(0)n(4 a) (0) 

4 4 ) ( )4 2) (±i) 

(with j = 0, . . . , n — 1; a = 1, 2 and a = 3 — a) 
4 1 i(0)Ai 4) (0)AW (0)A( 5 ) (±1) A( 3 ) (±2) 
41 (0)Af (0)Af (O)A(f) (±1) AW (±2) 
A«(0)4 5 )(0) 
Af(0)4f)(0) 

4 6 )(o)4 3 )(o)4 4) (o) 

4 5) (o) 



(B.16) 



Case D, (p=4n+6, n > 1): 



(n+1,2) (n+1,1) 




Bj B 3 K 2 Ki 

(0,2) (0,1) 



e B M = Rm Ba cosh 9- Yl $B a ,c*L c (6) (a = 1,2) 

ce{Bx,B2,Kt,K 2 } 
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ef\e) = 5 j0 [Rm Ka cosh 6- £ $ Ka4 * L d (9) + * Lf (9) 

\ de{B 1 ,B 2 ,K 1 ,K 2 } 13=1 

2 / n+1 \ 

- E *Stf * L f(°) - E ^" +2l ^ 2] * Lf{6) (B.17) 

/3=1 V fc=0 / 

with j = 0, . . . , n + 1, a = 1, 2. 
The Y-system is 

y Bl (±(p-2)) = A fi3 (0) 

y ft (±(p-2)) = A Bl (0)Ai 2 +1 (0)nA?Wn-4z)) 



i=0 



n+1 ^ jlDn+2] 



i=0 



(with j = 0, . . . , n + 1; a = 1, 2 and a = 3 — a) (B.18) 
B.3 Exceptional <f) 12 Y-systems 

The M PtP +i + 0i2 Y-systems are exceptional for p = 3, 4, 5 and 6. For p = 3, 4 and 6 
they are related to the exceptional Lie algebras E r with r = 8, 7 and 6. These systems 
are most conveniently written by departing from the conventions used elsewhere in this 
paper, and simply labelling the pseudoenegies as £i(0), . . . e r (9). The Y-systems are 



then gg]: 

y 3 (±(p - 2)) = f[ (A fc (0)) /fc/ U = 1, • • • ,r) (B.19) 

fc=i 

where is the incidence matrix of -Eg, -Ey or Eq for p = 3, 4 or 6 respectively. 

The Y-system for M^fi + 0i2 is new, and follows from the TBA equations given in 
§ra above. It is: 



^(±3) 


= A fi3 (0) 


3>B 8 (±3) 


= A Bl (0)A B5 (0) 


y BB (±3) 


= A fi3 (0) A K2 (±2)A Kl (0)A« (±1) A ( 2 ) (0) 


^B 2 (±3) 


= A^ 1 (±2)A( 1 )(±1)A^ 2 (0)A( 2 )(0) 


3V 2 (±1) 


= A flB (0)£W(0) 


3W±l) 


= A fl2 (0)£«(0) 


yW(±l) 


= £( 2 )(0)£^ 2 (0)£^(0) 


3^(±1) 


= £«(0) 



(B.20) 
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